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SUCTION  I 

rnnoDDcnoN/smisKi 


Recent  experiments  have  been  sponsored  by  the  Air  Force  Arnsnent  Labora¬ 
tory  (AFATL)  tn  measure  pressure  against  buried  structures  froa  high  length- 
to-dianeter  buried  explosives.  Hydrocode  calculations  were  required  to  help 
understand  anomalous  high  pressure  results  of  models  built  from  low  pressure 
measurements.  This  report  documents  the  several  calculations  performed. 

All  calculations  were  performed  using  a  model  of  Eglin  AFB  dry  sand  (40- 
percent  porosity).  The  calculations  predict  reflected  pressure  for  several 
incident  pressures,  free-field  pressure  in  the  sand  as  a  function  of  explosive 
size  and  reflected  pressure  on  a  plate  placed  in  the  soil  at  a  distance  of  6 
charge  radii. 

The  calculations  show  that  the  ratio  of  reflected  pressure  to  incident 
pressure  in  sand  is  a  function  of  the  incident  pressure,  reaching  an  almost 
steady  value  of  3.8  at  an  incident  pressure  level  of  30,000  PSI.  The  pressure 
ratio  is  slightly  less  than  3  at  an  incident  pressure  of  10,000  PSI. 


SECTION  II 


UDEOCODB  SMB  MODBL 


Calculations  were  performed  nsing  the  AFATL  version  of  the  TOODT  IV 
Lagrangian  wave  propagation  code  (Reference  1)  and  the  AFATL  BULL  Enlerian 
code  (Reference  2).  Both  employed  the  sane  sand  nodel.  The  nodel  is  based  on 
high  pressure  data  from  Van  Thiel's  Hugoniot  data  handbook  (Reference  3)  and 
low  pressure  data  taken  from  low  strain  rate  hydrostatic  neasurements  (Refer¬ 
ence  4).  The  sand  is  assumed  to  be  a  dry  quartz  sand  with  40  percent  poro¬ 
sity.  Its  initial  density  is  1.4  gm/cc  (100  lb/ft1). 

The  complete  pressure  versus  excess  compression  curve  is  seen  in  Figure 
1.  Arrows  indicate  whether  paths  can  be  used  for  loading  or  unloading  or 
both.  Yield  strength  is  assumed  to  be  zero  for  the  sand  (it  is  believed  to 
actually  reach  only  a  few  hundred  PSI  under  very  high  pressure  levels  and 
hence  is  neglected).  The  Gruneisen  ratio  for  sand  is  very  low  (perhaps  0.2  to 
0.3),  and  the  internal  energy  contribution  to  pressure  is  ignored. 


SECTION  III 
CALCULATIONS 


The  first  series  of  calculations  were  perforated  using  TOODY  and  square 
wave  pressure  inputs  at  various  levels  from  5,000  to  100,000  PSI.  The  calcu¬ 
lation  modeled  sand  backed  with  2— cm  thick  steel.  The  pressure  in  the  sand 
was  calculated  next  to  the  steel.  The  calculations  were  performed  in  one¬ 
dimensional  plane  geometry  using  symmetry  boundaries  in  the  two-dimensional 
TOODY  code. 

Figures  2  through  7  present  pressure,  pressure  impulse,  density,  and 
velocity  versus  time  for  incident  pressures  of  5,000,  10,000,  20,000,  30,000, 
50,000,  and  100,000  PSI.  Plotted  pressures  are  in  terms  of  kilobars.  These 
values  can  be  multiplied  by  14,500  to  convert  them  to  PSI.  For  example,  10 
kilobars  is  145,000  PSI.  Impulse  is  in  units  of  dyne-sec/cm*  (10»  dynes/cm* 
is  1  kilobar).  Drops  in  pressure  levels  in  the  plots  from  the  peak  reflected 
pressures  are  due  to  relief  waves  from  the  free  surface  of  the  steel  plate. 
The  plots  provide  Table  1  which  relates  the  reflected  to  incident  pressure 
ratio  to  incident  pressure. 

TABLE  1.  REFLECTED  TO  INCIDENT  PRESSURE  RATIO 


PjtPSI) 

PR(PSI) 

P  /P 

A  I 

5,000 

11,500 

2.32 

10,000 

28,320 

2.83 

20,000 

55,940 

3.30 

30,000 

105,275 

3.54 

50,000 

189,710 

3.79 

100,000 

391,450 

3.91 

These  ratios  are  plotted  in  Figure  8  and  demonstrate  that  PR/Pj  reaches  an 
almost  steady  value  of  3.8  at  the  50,000  PSI  incident  level. 
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To  apply  these  reflection  factors  to  the  experiment  of  interest*  one  mnst 
know  the  value  of  incident  pressure  level  in  the  sand  as  a  function  of  the 
size  of  the  explosive.  TOODY  calculations  were  performed  to  obtain  this 
information.  The  calculations  assumed  an  infinitely  long  cylinder  of  Composi— 
tion B  explosive.  Two  cylindrical  sizes  were  run  (4.0-cm  and  5.6-cm  radii). 
Figures  9  through  18  are  pressure  and  impulse  time  history  plots  from  the  5.6- 
cm  radius  calculation.  The  R  value  given  on  each  plot  is  the  initial  distance 
of  the  time  history  point  from  the  center  of  the  charge.  The  time  history 
points  move  with  the  flow.  Note  that  stations  1  and  2  are  inside  the  explo¬ 
sive.  Stations  in  the  sand  sometimes  exhibit  unrealistic  peak  pressure 
values,  especially  near  the  explosive.  These  large  peaks  could  be  damped  out 
with  additional  viscosity.  The  average  value  of  the  first  few  oscillations 
represents  the  true  pressure  at  the  station.  In  Table  2,  the  pressure  data 
has  been  adjusted  to  take  these  oscillations  into  account  by  extrapolating 
back  from  the  smooth  slope  of  the  pressure  curve  to  the  time  of  wave  arrival. 


TABLE  2.  FREE-FIELD  PRESSURE  LEVELS 

Explosive 

Station 

Peak 

BMl.bJL  ,(9») 

Location  (cm) 

Pressure  (kb) 

£SI 

5.6 

10 

300 

4,347.820 

14 

18 

260,870 

16 

14 

202,900 

20 

9.5 

137,680 

24 

5.3 

76,810 

28 

3.7 

53,620 

32 

2.7 

39,130 

The  calculation  with  the  explosive  radius  set  to  4  centimeters  yields  the  same 
incident  pressure  levels  when  plotted  versus  R/Rgg  where  R  is  the  distance  to 
the  station  measured  from  the  axis  of  the  explosive  and  Rgg  is  the  explosive 
radius.  The  scaling  factor  for  an  infinitely  long  explosive  then  is  R/Rgg. 


Figure  19  presents  incident  pressure  versus  R/Rgg  data  from  these  two  calcula¬ 
tions.  The  pressure  from  R/Rgg=5  to  7  (and  it  is  assumed  to  extend  to  8)  can 
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Figure  13.  Station  5  Free-Field  Pressure  and  Impulse 


m* 


be  approximately  fit  with  the  equation: 

-0.376 

Pj  -  23 .4  e 


(kb) 


This  equation  is  not  valid  for  R/Rg£  values  less  than  5.  For  normal  incidence 
of  the  sand  shock  vave  with  a  thick  steel  plate,  the  reflected  pressures  given 
in  Table  3  can  be  expected  at  scaled  distances  of  5  to  8.  The  pressures  are 
calculated  from  using  the  incident  pressure  equation  above  with  Pg/Pj  values 
from  Table  1. 


TABLE  3.  REFLECTED  PRESSURE  VERSUS  SCALED  DISTANCE 


r/rbe 

5 

6 

7 

8 


PjtPSI) 

51.750 
35,530 
24.395 

16.750 


PR(PS1) 

196,650 

127,910 

83,080 

52,720 


The  next  calculation  utilized  the  BULL  code  and  was  designed  to  investi¬ 
gate  the  variation  of  pressure  along  a  flat  plate  established  at  a  scaled 
distance  of  6  from  the  center  of  a  5.6-cm  radius  Composition  B  charge.  Figure 
20  is  a  density  contour  plot  showing  the  undetonated  explosive  in  sand.  The 
two  axes  (7  vertioal  and  R  horizontal)  are  planes  of  symmetry  and  the  equa¬ 
tions  solved  are  those  applying  to  plane  geometry.  The  charge  is  therefore  an 
infinitely  long  bare  explosive.  The  wall  is  a  rigid,  perfectly  reflecting 
plane  at  an  R  value  of  33.6  cm  (R/R  *6.0).  It  does  not  show  up  in  the 

nil 

contour  plots  until  impacted  by  a  shock  wave  at  approximately  200  micro¬ 
seconds.  Figures  21  through  47  are  density  contours,  pressure  contours,  and 
pressure  histograms  along  the  Y  (vertical)  axis  at  approximately  50-micro¬ 
second  intervals  to  400  microseconds.  Note  that  the  pressure  histogram  des¬ 
cribes  free-field  pressure  since  it  presents  values  along  the  vertical  axis. 
Reflected  pressures  will  be  seen  in  time  history  plots  presented  later  in  this 
report.  The  discontinuities  in  pressure  halfway  up  the  pressure  histograms 
are  oscillations  at  the  explosive/sand  interface.  Numbers  on  the  right  side 
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Figure  20.  Initial  Setup:  Composition  B  in  Sand 
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Figure  21.  Density  Contours  at  50  Microseconds 


PRESSURE-  DYNES/01-30  IX IB  *J 

OTI  HULL  CALC  -  5.6  CM  RHE  -  RIGID  WALl  AT  33.6  CM  FROM  ORIGIN 
TIME  50.000  U5EC  CYCLE  334.  PROBLEM  2127. 1525 


Figure  23.  Pressure  Histogram  at  50  Microseconds 
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Figure  24.  Density  Contours  at  100  Microseconds 
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Figure  26.  Pressure  Histogram  at  100  Microseconds 
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Figure  28. 


Pressure  Contours  at  150  Microseconds 
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Figure  30.  Density  Contours  at  200  Microseconds 
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Figure  31.  Pressure  Contours  at  200  Microseconds 
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Figure  32.  Pressure  Histogram  at  200  Microseconds 
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Figure  34.  Pressure  Contours  at  225  Microseconds 
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Figure  37.  Pressure  Contours  at  250  Microseconds 
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Figure  38.  Pressure  Histogram  at  250  Microseconds 
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Figure  39.  Density  Contours  at  300  Microseconds 
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Figure  40.  Pressure  Contours  at  300  Microseconds 
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Figure  41.  Pressure  Histogram  at  300  Microseconds 
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Figure  43.  Pressure  Contours  at  350  Microseconds 
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Figure  45.  Density  Contours  at  400  Microseconds 
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Figure  47.  Pressure  Histogram  at  400  Microseconds 


and  across  the  tops  of  sone  of  these  figures  refer  to  zone  numbers  in  the 
calculation.  The  (1*1)  on  the  histograms  means  that  the  data  apply  to  the 
first  radial  zone.  This  report  has  used  R  and  ¥  as  axis  descriptors.  Minimum 
and  maximum  contour  values  in  the  figures  are  indicated  as  existing  at  X  and  ¥ 
positions.  In  this  case,  the  X  axis  is  the  same  as  the  R  axis. 

The  first  shock  wave  is  seen  to  hit  the  wall  between  200  and  225  micro¬ 
seconds  (see  Figures  29  and  32).  Figure  34  shows  that  a  mild  Mach-stem  has 
developed  along  the  wall  at  the  intersection  of  the  reflected  and  incident 
waves.  This  behavior  continues  as  seen  in  subsequent  pressure  contour  plots, 
although  one  must  mentally  draw  in  the  incident  pressure  curve  from  the  lead 
location  on  the  pressure  histograms,  since  contour  values  chosen  by  the  plot¬ 
ting  program  do  not  include  it  after  225  microseconds.  The  establishment  of  a 
Mach-stem  leads  to  slightly  higher  pressures  at  positions  away  from  the  point 
closest  to  the  explosive.  This  effect  is  seen  in  the  pressure  and  impulse 
time  histories  recorded  at  locations  near  the  wall. 

Fignres  48  through  55  are  pressure  and  pressure  impulse  time  histories 
for  stations  located  in  the  sand  next  to  the  rigid  wall.  The  stations  are 
fixed  with  respect  to  the  wall  and  do  not  move  as  the  sand  moves.  For  sta¬ 
tions  beyond  the  first  two  (48  and  49),  judgement  is  again  required  to  deter¬ 
mine  the  actual  peak  as  distinguished  from  the  initial  oscillations. 
Pressure-time  slopes  from  the  smooth  part  of  each  curve  were  extrapolated  back 
to  the  time  of  arrival  of  the  wave  to  generate  Table  4. 

Free-field  pressures  at  various  distances  from  the  explosive  can  be 
determined  by  smoothing  the  peak  pressures  in  the  pressure  histogram  curves  up 
to  300  mioroseconds.  Beyond  this  time  the  reflected  pressure  from  the  wall 
has  reached  the  Y  axis  and  disturbed  these  values.  Examination  of  the  histo¬ 
gram  values  shows  excellent  agreement  with  the  free-field  data  in  Table  3. 
The  incident  peak  pressure  then  at  the  wall  is  35,530  PSI.  The  reflected 
pressure  at  the  point  of  normal  incidence  is  158,116  PSI,  yielding  a  PR/Pj 
ratio  of  4.45.  It  is  believed  that  this  reflected  pressure  ratio  is  larger 
than  that  obtained  previously  because  of  the  nature  of  the  absolutely  rigid 
wall  in  this  oaleulation.  The  steel  in  the  previous  calculations  was  able  to 
move  slightly  and  relieve  peak  pressure. 


53 


6  CM  RTC  -  RIG10  WALL  AT  33  6  CM  rROM  ORIGIN 


STA  |  R«  33. >.500  v.  .mOO  MIN-  I  .  39«E  03  MAX-  I.091E-01  LAST-  J.373E-0I 


- - - - ,  -  ,  - - —x/-  r 

4.000E-0I  B.OOOC-OI  1.500E-0?  I.600E-OS  S  OOOE-O?  S.400E-02  £  90 OE *05  3. 

T  IME(  USEC) 

OTI  HEJLL  CALC  -  3.6  CM  RTC  -  RIOIO  WALL  AT  33.6  CM  FROM  ORIGIN 


STA  1  R-  33.4500  Y-  .1500  MIN*  I  .081E-00  MAX*  4.H8E*05  LAST-  4.II8E*05 


r  r  me  c  usec ) 


Figure  49.  Station  2  Pressure/Impulse  at  Rigid  Wall 


T MLB  4.  REFLECTED  PRESSURE  AT  THE  TALL 

STATION 
NUMBER 

1 

2 

3 

4 

5 

6 

7 

8 

Examination  of  Table  4  data  shows  that  reflected  pressure  increases  away 
from  the  normal  incidence  point  and  then  decreases  again.  This  is  dne  to  the 
effect  of  the  Mach-stem  formation.  Table  5  presents  incident  pressure  com¬ 
puted  from  the  previous  scaled  range  equation  using  distance  from  charge 
center  to  the  station.  The  incident  pressure  is  seen  to  drop  from  35,530  PSI 
at  atation  1  to  13,045  PSI  at  station  8.  Using  this  incident  pressure  value 
and  actual  station  pressure,  values  of  Pg/Pj  were  calculated.  Multiplying 
this  ratio  by  the  cosine  of  the  angle  the  flow  makes  with  the  normal  to  the 
wall  at  the  station  location  further  smooths  the  data. 


TABLE  5.  REFLECTED  PRESSURE  VS  ANKLE 


STATION 

Y  POSITION 

ANGLE 

Jin 

h!h 

rR 

NUMBER 

—(.cm _ 

(PEG) 

PI 

1 

0.15 

0 

35,530 

4.45 

4.45 

2 

5.15 

8.75 

34,953 

4.88 

4.82 

3 

10.15 

16.88 

32,438 

5.36 

5.13 

4 

15.15 

24.37 

28,813 

5.98 

5.45 

5 

20.15 

31.06 

24,642 

6.47 

5.54 

6 

25.15 

36.94 

20,419 

7.46 

5.96 

7 

30.15 

42.02 

16,490 

8.79 

6.53 

8 

35.15 

46.42 

13,045 

5.34 

3.68 
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The  calculations  presented  in  the  previons  section  nnst  be  applied  to 
test  results  with  caution  since  they  all  assuae  infinitely  long  cylindrical 
charges.  End  effects  will  obviously  have  sons  influence  if  the  charge  is  of 
finite  length.  However,  it  is  believed  that  they  can  be  used  to  help  inter¬ 
pret  test  data.  The  calculations  clearly  show  that  very  large  reflections  can 
be  expected  at  scaled  distances  of  6  or  less  froa  the  charge  and  that  pressure 
will  vary  across  a  steel  plate  placed  in  the  flow  due  to  varying  distance  from 
the  charge,  varying  flow  angle  with  respect  to  the  plate,  and  Mach-stea 
formation. 
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